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We study the effect of mass imbalance on the phase diagram of a two-component fermionic mixture 
with attractive interactions in optical lattices. Using static and dynamical mean-field theories, we 
show that the pure superfluid phase is stable for all couplings when the mass imbalance is smaller 
than a limiting value. For larger imbalance, phase separation between a superfluid and a charge- 
density wave takes place when the coupling exceeds a critical strength. The harmonic trap induces 
a spatial segregation of the two phases, with a rapid variation of the density at the boundary. 
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The remarkable advances in handling ultra-cold atomic 
gases have given birth to the new field of "condensed mat- 
ter physics with light and atoms" . Cold atoms in optical 
lattices, with tunable and controllable parameters, have 
been studied in many different contexts (for reviews, see 
Mixtures of two-component atoms with different 
masses (e.g ^Li, "^"K) introduce an additional parameter, 
namely the difference between the hopping amplitudes 
associated with each species in the optical lattice. This 
may affect the stability of the possible quantum phases 
or even induce new ones. Recently, a phase diagram has 
been worked out in the one-dimensional (ID) case 0. 

In this article, we consider such fermionic mixtures in 
higher dimensions, with an attractive on-site coupling. 
Using analytical and numerical techniques, we establish 
a ground-state phase diagram as a function of coupling 
strength and mass imbalance, in all regimes of couplings. 
We also consider the experimentally relevant effect of the 
trap potential, which is shown to induce a spatial segre- 
gation between superfluid and density-wave phases. 

Under conditions discussed, e.g., in Refs. P, H, 01, 
fermionic mixtures are described by a Hubbard model: 

H = - J2 + h-c-) - |U| ^ ni^nu (1) 

The (pseudo-) spin index a refers to the two different 
species. Feshbach resonances between ^Li and are 
currently under investigation 0], and would allow for 
an attractive interaction with a tunable strength, as as- 
sumed in ID). For an example of hetero-atomic reso- 
nances in the boson- fermion case, see e.g [gj. In the 
following, a bipartite optical lattice made of two inter- 
penetrating {A, B) sublattices (such as a cubic lattice) 
is considered. For simplicity, we consider an equal num- 
ber of atoms for each species, leaving for future work the 
study of imbalanced populations. 

In order to study the ground-state phase diagram of 
model (H)), we use dynamical mean field theory at zero 



temperature (DMFT) [7], together with analytical mean- 
field calculations for both weak and strong coupling. Let 
us anticipate the DMFT phase diagram of the uniform 
system, displayed in Fig. [TJ When the fermions have the 
same mass, the ground-state is a superfluid (SF) for all 
J7| . A competing ordering exists, namely a charge den- 
sity wave (CDW), considered here in the simplest (com- 
mensurate) case in which the charge is modulated with 
an alternating pattern on the A and B sublattices. At 
half-filling {{n^ -\- ni) = 1), it is well known that the SF 
and CDW states are degenerate. This no longer applies 
in the 'doped' system away from half-filling: for equal 
masses, the SF phase is stabilized by doping for all \U\ , 
but a large mass imbalance favors the CDW phase over 
a SF state in which the Cooper pairs must be formed by 
fermions with different mobilities. Hence the SF/CDW 
competition becomes more interesting in the presence of 
mass imbalance. As displayed on Fig.[Tl we find that the 
uniform system has a SF ground-state for all values of 
\U\ as long as the mass imbalance z = (t| — t|)/(tf -I- 1|) 
is smaller than a limiting value Zc (which depends on the 
average density). For z > Zc, a (first-order) phase bound- 
ary is crossed as |J7|is increased, beyond which the un- 
form system undergoes a phase separation (PS) between 
a SF and a CDW phase. As discussed later in this paper, 
this implies that, in the presence of a harmonic trap, the 
CDW and SF phases may both exist in different regions 
of the trap. 

DMFT is a quantum generalization of classical mean- 
field theories, which takes the full local quantum dy- 
namics into account, while spatial fluctuations are ne- 
glected. It maps a lattice model onto an effective 'quan- 
tum impurity model' (a single interacting site which hy- 
bridizes with an uncorrelated bath), subject to a self- 
consistency condition Q. DMFT and its extensions 
have been used to study the attractive Hubbard model 
with equal masses Q. It is convenient to work with 
Nambu's spinors — {c^ , C| ) . The key quantity con- 
sidered in DMFT is the local (on-site) Green's function: 
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^i''') = (Tr'0i(T)f/'i-'~(O)) and its Fourier transform for 
imaginary frequencies: 
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The superfluid order parameter is then given by Asf = 
(c,|C,x) = F(r = 0) = J2^F{iuj). In the CDW state, 
the local Green's function takes different values {Ga and 
Gb) on each sublattice. The CDW order parameter is 
the difference of densities on each sublattice: Acdw — 
{uA ^ tlb)- The self-consistency conditions of DMFT 
relate the (frequency-dependent) 'Weiss fields' 
entering the effective 'impurity model' on one sublattice, 
to the Green's functions, through 7]: 

= + P'-TGB{A){i(^)T, (3) 

in which T = diag[t|, — tj and fi = diag[/Z|, —iJ-i] are di- 
agonal matrices associated with the hopping and chemi- 
cal potential of each species. As written, ([3]) assumes for 
simplicity a semi-circular density of states, but is easily 
generalized to an arbitrary lattice. Eq.® allows for the 
study of both SF and CDW orders, and for their possi- 
ble coexistence. The ground-state energy of the different 
phases is evaluated as (H) — {K) + UJ2i{''T'i1^il)j '^it^^ 
the kinetic energy (K) in the SF and the CDW phases 
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reading, respectively: {K)sf = 
F^iLo)] and (iC)cDW = 

We performed DMFT calculations [l^ spanning the 
whole range of coupling \U\ and imbalance z. We fo- 
cused on the vicinity of half-filling, and found the phase 
diagram of the uniform system (Fig. [T]) to be qualita- 
tively independent of the 'doping level' S — {n^ + — 1). 
For small enough values z < Zc{S) of the mass imbalance, 
a pure SF solution is stable for all |C/|. In contrast, for 
z > Zc, the pure SF phase is stable only for small in- 
teractions (below the line drawn in FiglT]). Above this 
line (which depends on i5) , the pure SF solution becomes 
unstable towards phase separation between a SF and a 
CDW phase. (Note that we did not find a homogeneous 
CDW solution out of half-filling, except at z = 1). This 
means that it is more convenient to separate the system 
into a fraction 1 — x with CDW order and i5 = 0, and a 
fraction x with SF order accommodating the rest of the 
particles. This conclusion is reached by minimizing over 
X the expression £'ps(a;) = (1 — x)£'cdw + xEsf- We 
note that the SF phase is more stable than in the ID 
case (in which nesting favors a CDW with Q — 2fci?). 

In the following, we describe analytical mean-field cal- 
culations for both weak and strong coupling which help 
in understanding the DMFT phase diagram established 
numerically. We first present a strong- coupling analysis, 
which holds for \U\ ^ t|,t|. In order to analyze this 
limit, we find it useful to resort to a particle-hole trans- 
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FIG. 1: Phase diagram of the uniform system in the {z, \ U\) 
plane, obtained from DMFT. Below the curves (displayed here 
for two 'doping' levels 5 = n — 1 = 0.05, 0.1), the superfluid is 
stable. Above the curves, the system is phase-separated into 
a half-filled CDW and a SF. The arrows indicate the strong- 
coupling values obtained analytically. The dotted lines are 
the weak-coupling mean- field approximation (see text). \U\ is 
normalized to the bandwidth W of (skt + eki)/2. 



formation (Table [T| and work in the repulsive-f7 frame- 
work. We emphasize that we are not switching to truly 
repulsive interactions, but we simply exploit a mapping. 
Under this mapping, our model is transformed, at large 
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TABLE I: Particle-hole transformation mapping the U < Q 
model with (n|) = {n\) onto a half-filled U > Q model with a 
magnetic field. 

\U\ ^ t|, t|, into an XXZ quantum spin- 1/2 model Hi]: 
H = jY^ Si.Sj +iJY^ SfS'j - h'^{2St - to) (4) 

in which S = ^d^dafjdp, J = 4f|t|/|[/| and 7 ~ 
(t| — /2ttti = 22:^/(1 — z^). Hence, the mass im- 
balance turns into a spin exchange anisotropy. The uni- 
form magnetic field h corresponds to the original chem- 
ical potential /i — |C^|/2 and the magnetization to the 
doping 5 (cf. Table. |T|. The mean- field approach 
amounts to treat the spin variables as classical, and min- 
imize the energy over the angles 0a, Ob describing the 
orientation of the spins in the two sublattices. The en- 
ergy per site reads (with C the lattice connectivity and 
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CA,B = cos9a,b,sa,b =sin9A,B)- 

= SASB+'^J{^+l)cACB~^[cA + CB-'2rn] (5) 

The phase diagram is characterized by the competition 
between the xy spin-density wave (SDW^;!,) with order 
parameter l^^y = {{—^ySf) (corresponding to SF order- 
ing for U < 0), and Neel order (SDW^), = ((-l)'S'f ) 
(corresponding to CDW). The solution changes accord- 
ing to the magnetization m of the system (i.e. the doping 
of our physical model). The m vs. h curve has a discon- 
tinuity of amplitude rric = ^7/(7 + 2) = z. For m — 
(half-fiUing 6 = 0), a SDW^ (CDW) state is obtained. 
For m G [rric, 1], the homogeneous SDWa;^ (SF) state is 
stable, while for < to < TOc phase separation takes 
place between the two types of ordering. Thus, when 
working at fixed magnetization (corresponding to fixed 
doping S), one finds a SF for z < Zc ^ rn = S and phase 
separation for z > Zc = S. This strong coupling value 
(indicated by arrows on Fig. [T]) agrees very well with our 
DMFT results. 

We now turn to the opposite weak-coupling limit. We 
decouple the interaction term in the SF and the CDW 
channels, and determine the regions of stability of each 
phase. We first consider the BCS decoupling of the 
interaction, introducing the order parameter Abcs = 
{\U\/N) J^ki'^ki^ki) ^'^ make the Hamiltonian quadratic. 
In Nambu formalism it reads: 



^^Bcs = X! 
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Here, jl^ = - Un^a, ^ko- = £ko- - Mo- and Eq = 
EkCki + N\U\n'^ni + iVA|cs/|[/|. The diagonaliza- 
tion of (O yields the Bogoliubov modes with eigenvalues 
Et = ±(CkT 



■ai)/2+V(aT+ai)V4 + A|cs- 



Defin- 



ing new variables = (Ckt + Cki)/2, ji ^ {jl^ + Ai)/2, 
the usual form of the BCS gap equation is recovered, and 
tells us that the normal state is always unstable toward 
SF ordering. Analogously, we can decouple the interac- 
tion in the CDW channel defined by the order parameter 
A. = (|C/|/iV)Ek(4+Q.Ck.) with Q = (^,... ,7r). In- 



troducing the spinor -0] 
Hamiltonian reads 
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the mean-field 



CDW 



keRBZ.o- 



Ekcr — Mcr —Act 
— Act — Ekcr — Act 



V^kCT + Eq 

(7) 



with Eq = 7VA|A|/|;7| + N\U\n^ni. It is readily diag- 
onalized, with eigenvalues: E^^ — ±-\/e^~rAj — jla- 
This yields two self-consistency conditions: 
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FIG. 2: Phase diagram for S — 0.05 from weak-coupling mean- 
field (whose validity is questionable above the dotted line) -see 
text-. For simplicity, a square density of states was used here. 



At a fixed value of the chemical potential, these CDW 
equations have the following solutions: i) for all \U\ and 
z, a normal solution with Acdw = 0, which is unstable 
towards SF ii) for large enough |J7|, a half-filled (commen- 
surate) CDW and iii) for large values of z, close to 1 a 
homogeneous CDW solution is also found with a density 
different from unity {S =/= 0). 

We first compare the ground-state energies of two 
mean-field solutions: the homogeneous SF, and the 
SF/half-fiUcd CDW phase-separated solution obtained 
from a Maxwell construction. The resulting phase 
boundary (Fig. [T]) is seen to be qualitatively reasonable, 
and even quantitatively accurate (in comparison to the 
numerical DMFT result) for some intermediate range of 
z. Indeed, the weak-coupling mean-field is justified when 
\U\ < t|,tx, i.e. \U\/W < (1 - z). In Fig. El we perform 
a more detailed comparison of the ground-state energies 
of three mean-field solutions: the homogeneous SF, the 
phase separated SF/CDW, and the homogeneous CDW 
with S ^ (when it exists). This comparison yields a 
small region of parameters, for large z, in which a ho- 
mogeneous CDW with a density different from one atom 
per site is stable. This pure CDW pocket might be unsta- 
ble to SF ordering, yielding a candidate for a supersolid 
region, but we did not check this explicitly. Anyhow, 
this solution is stabilized in a region where the reliability 
of weak-coupling mean-field is questionable. The lack of 
such a solution in DMFT may lead to the conclusion that 
the CDW (or supersolid) is an artefact of weak-coupling 
mean-field, but it must also be noted that the DMFT 
solution has a finite numerical resolution, and that the 
energetic balance involved is very delicate. Hence, we 
cannot reach a definitive conclusion on this issue. Ex- 
actly for z = I the "down spin" atoms are no longer 
mobile and we have a Falicov-Kimball model, which has 
a pure CDW ground state ■ 

We finally discuss the effect of the trap potential. For 
simplicity, we perform an explicit calculation only in the 
strong coupling limit, using again the particle- hole trans- 
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formation (Table [T| and considering the effective spin 
model ([4]). A harmonic trap potential yields a position- 
dependent chemical potential which corresponds, under 
the particle-hole transformation to a spatially varying 
magnetic field h{r) — h — h^r^/R^. Here i?o is the radius 
of the circular trap, /iq = and h = ^ — \U\/2 is 

related to the chemical potential at the center of the trap, 
which must be adjusted so that the local density n{r) in- 
tegrates to the total number of atoms. We start from 
a local density approximation (LDA), and also compare 
with a Monte Carlo solution of the strong-coupling model 
in the presence of h{r). As described above, the strong 
coupling analysis of the uniform system yields a critical 
magnetic field (chemical potential) at which m(/i^ is dis- 
continuous. For \h\ < he — JCV^h + 2) = T^'fuY' '^^ 
have a SDWz (CDW) phase, otherwise we have a SDWa;^ 
(SF) phase. Within the LDA approximation, this implies 
that in a region where \h{r)\ is smaller (resp. larger) 
than he we locally observe SDWz/CDW ordering (resp. 
SDWa;y/SF). According to the values of the parameters 
h and ho, and noting that h — ho < h{r) < h, one finds 
several different regimes: 

(i) h—ho > he or h < —he- The trap potential is always 
larger than he, or smaller than —he, so that the system is 
in a ST)Wxy (SF) phase everywhere inside the trap, and 
the density profile varies smoothly, (ii) h > he and \h — 
ho\ < he', in this case, h{r) > he inside a circle of radius 
i?i = Ro\/ {h — he)/ho centered at r = 0. Hence, one has 
phase separation into two distinct regions: SDWa;j,(SF) 
ordering within this circle, and SDW^ (CDW) in the outer 
ring (Fig. [3l left panel), (iii) h ~ ho < —he and \h\ < he: 
we find again phase separation, with the opposite spatial 
arrangement. The SDWj;j,(SF) part is stable out or a 
circle of radius i?2 = + he) /ho, inside which there 

is a SDW^(CDW) phase (Fig. [3 middle panel), (iv) h > 
he and h — ho < -he- Then, the magnetic field profile 
crosses both he and -he, so that there are three spatial 
regions: R < Ri where we find SDWj;y(SF), then the 
ring i?i < r < i?2, where SDW2(CDW) establishes, and 
finally an outer ring r > R2 with SDWa;j,(SF) ordering 
(Fig. [3l right panel). 

In the three last cases (ii-iv), in which phase separation 
occurs, the LDA approximation predicts a jump of the 
magnetization at the phase boundaries i?i and R2, corre- 
sponding to a jump of the density in the original U < 
model (see also [111). In order to test this prediction 
and assess the validity of LDA, we performed a classical 
Monte Carlo simulation of model (|4|) in the presence of 
a spatially dependent field h{r). For simplicity, this test 
was performed in a one- dimensional geometry. We find a 
remarkable agreement between the LDA density profiles 
and the Monte Carlo solution, which confirms that very 
sharp variations of the local density indeed takes place 
at the boundary between domains in cases (ii-iv). 

In conclusion, for attractive interactions, the domi- 
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FIG. 3: Density profiles and domains with different orderings 
inside the trap (bottom panels), as discussed in text. The top 
panels show how the trap potential intersects the character- 
istic values of the chemical potential in each case. 



nant effect of the mass imbalance is to induce a compe- 
tition between superfluid and density-wave ordering. In 
the presence of an inhomogeneous trap potential, both 
phases can be stabilized in different regions of the trap, 
with rapid variations of the local density at the phase 
boundaries. We note finally that, in the case of the 
^Li/^^K mixture, a simple estimate shows that the mass 
imbalance z can be varied over a large range by chang- 
ing the lattice depth Vo/Eji (z ^ I at small Vo/En and 
z ~ 0.9 for Vo/Eji ~ 15), so that the effects discussed in 
this work may indeed be observable in this system. 
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